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Abstract 

We consider the problem of tunneling escape of particles from a multiparticle system confined 
within a potential trap. The process is nonlinear due to the interparticle interaction. Using the 
hydrodynamic representation for the quantum equations of the multiparticle system we find the 
tunneling rate and time evolutions of the number of trapped particles for different nonlinearity 
values. 
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The recent advances in the experiments on real Bose - Einstein condensates {BECy^^ 
and nonhnear optical waves^ have generated a huge body of works on the theoretical side. 
The standard approach on the theory side was the use of the nonlinear Schrodinger equation, 
with a potential term depending on the particle density (also known as the time dependent 
Gross-Pitaevskii equation^'^'^*^) . However, as it became clear early in the investigation, the 
time behavior of such equations was very complex and rich. The phenomena of coherence^^, 
macroscopic tunnelingii, vortex formation , ^^i^^i^^^i^^ instabilities, focusing and blowup are 
all new concepts which are related to the nonlinear nature of the systems. Most of the 

;ic problems are so far being dealt with by a 



analysis of these hard and fundamental analy 
combination of numerical schemes (e.g. Ref. ^) 



and finite dimensional phenomenological 



models. The quest for a theory that can adequately give the relevant nonlinear effects in 
the time dependent regime is therefore of major current interest. 

In this letter we analyze the problem of resonance and tunneling of N-body Quantum 
systems by solving the corresponding nonlinear problem in the leading relevant approxima- 
tion. Our approach combines ideas from Many body theory, nonlinear partial differential 
equations and resonance theory in Quantum Mechanics to offer a new unified approach to 
finding tunneling times for both linear and nonlinear systems. We focus on the problem of 
a BEG droplet in a potential well, that can tunnel through a finite barrier from the trap. 
We find the leading nonlinear corrections to the tunneling rate. 

A droplet of a large number of atoms, A^, with boson statistics is confined by an external 
potential Vext(R)- The multiparticle wave function describing the quantum state of such a 
system satisfies the Schrodinger equation 

Starting from this equation we may derive equations of motion for moments of the single 
particle Wigner function p\y{p,£, R, t).^-^ We concentrate here on the condensate behavior, 
which is supposed to be separated from the dissipative excitations (a detailed discussion of 



ih- 



this separation can be found in Ref. llTl ). The derivation can be made gauge invariant and 
its details are presented in the papersi^ where many relevant references can be also found. 
The most important of the moments are the density distribution of the particles 



and the velocity field 



for which the continuity equation 

dt 

and the Euler-type equation 



V(v(R,t)p(R,t)) = (2) 



^^^^^ + (v(R,t) ■ V)v(R,t) = -lv(K..(R) + V,^{R) + K//(R)). (3) 

are obtained. The velocity field is not necessarily potential and three dimensional cases 

^2 y^\/p(R~t) 

allows for vortices. Here the "quantum potential" V^^fR) = ' - accounts for 

the quantum character of the liquid in the droplet, whereas Ve//(R) is the contribution due 
to the inter-particle interaction. In many cases, the interaction potential Ve//(R) can be 
represented as a functional of the particle density p(R, t). For example, if we assume that 
the inter-particle interaction is point-like, Vint(Ri — Rj) = |A(5(Rj — Hj) then the direct 
calculations, whose many aspects are parallel to those presented in Ref.-^, result in the 
lowest order in A in Ve//(R) = Ap(R, t). Then equations and (jHI) become equivalent 
to the time dependent Gross-Pitaevskii equation. Higher order term in A (to be discussed 
separately) will produce a correction to the above effective potential Ve//(R) as well as a 
dissipative contribution due to the interaction with the excitations above the condensate. 

In this presentation we will limit ourselves to discussion of the problem of tunneling 
evaporation of the droplet kept within a one-dimensional potential depicted in figure ^ It 
means that we will have to consider the one dimensional versions of equations Q and ©, 
which read now 

d 

pt{x,t) + —[p{x,t)v{x,t)] = 0. (4) 



and 

1 d 

Vt{x, t) + v{x, t)Vx{x, ^) = 



Vecctix) == — VpOM) + Ueff{p{x,t)) 

^/p{x,t) 2m 



(5) 



In order to solve the problem we first introduce an auxiliary confining potential Vo(x), 
which coincides with the real potential Vext{x) for small x but differs from it (dashed line 
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FIG. 1: Potential well keeping the droplet. The dashed line shows the auxiliary potential Vo{x) at 
larger x. The dash dotted line corresponds to the energy E{N) for a given number of particles A^. 

in fig. at large x. Stationary states are possible in such potential with the density 
distribution following from equation (0) at v{x, t) = and with Vq{x) substituted for Vext{x). 

Vo{x)^p{x,t) - — — ^p{x,t) + Ueff{p{x,t))^/pix,t) = E{N)^p{x,t) (6) 

E{N) is the single particle eigen-energy in a droplet with particles. One readily 
recognizes a stationary Schrodinger-type equation for the single particle wave function 
4>{x,t) = a/ p{x, t) and with the effective potential Ueff{p{x,t)) accounting for the many 
body effects. In case of a contact interaction between the particles equation (jHl) becomes 
celebrated stationary Gross-Pitaevskii equation. 

The solution of equation © at a given initial number of particles Nq is considered as 
an initial state of a droplet in the real external potential Vext{x)- In this case the density 
distribution p{x,t) becomes time dependent and its dynamics is governed by equation 
Here we will follow the time evolution of the number of particles N{t) within the droplet, 
defined as N(t) = dxp{x,t) and with xq chosen as the exit point of the potential Vext{x) 
at a given energy E{N) (see fig. Q). N{t) at t = nearly coincides with the initial number 
of particle Nq since only exponentially small part of the particle density remains outside the 
potential well at x > Xq. 
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Equation © can be solved if we assume that the velocity field and the density distribu- 
tion change adiabatically slow with time. The applicability of this approximation will be 
discussed below. Then we may neglect the time derivative Vt{x,t) = in the rhs of ©. 
Using the fact that initially the density distribution p{x, t) satisfies equation (jH)) we may 



find the velocity field v{x,t) = ^ ^AV{x) where AV{x) = Vq{x) — Vext{x) 

Now we find the number of particles within the well Vext{x) by integrating the continuity 
equation from to the exit point xq. According to our approximation both the density 
and the velocity fields vary slowly with time only due to the time variation of the number 
of particles N{t). The exit point xq is connected with eigen-energy E{N) of the tunneling 
particles, so that the velocity becomes v{xo,t) = ^^[Vq — E{N{t))] and we get equation 

~ 1 

t = - dN =. (7) 

p{xo)J-[Vo-Em 

which determines implicitly the time dependence of the number of particle N{t). 

In order to estimate the particle density p(xo) at the exit point the particle density at 
the eigen energy E{N) is written as p{x; N) = Ng{x;N), with g{x; N) normalized to one. 
p{x, N) satisfies Schrodinger equation (jHI) but we need to know its value only at the exit 
point where it is exponentially small and we may neglect the inter-particle interaction. As 



a result we get ^/ g{xo; N{t)) oc e""^" with a = j:^/2m{Vo - E{N)). 

If the dependence of the energy E{N) on the number of particles in the droplet is known, 
equation ((Tj) can be solved and the time dependence N{t) of the number of particles can be 
found. As an example we carry out the calculation for the contact inter-particle interaction 
when equations (j^ and (0) correspond to the Gross-Pitaevskii equation. At not too large N 
we may approximately assume a linear dependence of the energy on the number of particles 
E{N) = £'(0) — XN— and represent the integral (|7j) in the form 

1 n Q^ooxoy^ 



n^ot = - dy (8) 

where = 2^, Ngp = ^. ao = y2m{Vo - E{0)) and Tq = ^g{xo; 0). 
Using the property 

1 ^^^^ y + ^/fTl-1 
y^l + y^ dy y+^Jy^ + 1 + 1 
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the integral (jS} is taken by parts and after that the slow logarithmic function is taken out 
of integral. As a result we get 



-Tot 
2 " 



y+^/y^ + 1-1 _ 1/0 + Vvi + 1 - 1 



g^aoxoy _ (9) 



y+,/y^ + 1 + 1 Z/0 + Vl/o + 1 + 1. 
This equation can be solved iteratively. First we assume that y = yo in the exponential 
function equation and get 

2f{t;yo) 



where f(t]y) = v{y)e a'""^ ^ ^ and viy) = ^Ltv^^L-^, or ?/ = ^'^l^L ■ At the next 
iteration we substitute solution (fTIHl into the function f(t; y) and repeat the procedure, 



2mA \l - r-(t; vi(t)) 

The function i/ in the above equations is always less than one. v tends to zero in the 
linear limit, when A 0, and one gets the exponential decay N{t) = Nqc'^^^. At a large 
nonlinearity parameter A, the function u may be close to one and the initial decay strongly 
deviates from the exponential behavior. However, u diminishes in the time course and at 
large time the asymptotical behavior of the number of particle N{t) becomes exponential. 

Before discussing the results it is worthwhile to estimate the validity of the adiabatic 
approximation made while solving equation (0). For this we need to compare the time 
derivative of the velocity Vt near the exit point xq, neglected in the above calculations, with 
the right hand side of equation The time dependence of the velocity results mainly from 
the time dependence of the exit point, so that 

vt = r^Kxt(a;)U=xo ~ T^T^o ^2 

mv[Xo) mv[Xo) I 

with I being the typical length characterizing the potential Vext{x). Since xq is the solution 
of equation Vext{x) = E(0) — \N(t) we get x ~ ^if"^- The right hand side of equation (jSj) 
can be estimated as ~ Vo/ml so that comparing it with the time derivative of the velocity 
field (fT^ we need to check the inequality 

vtml xo NXN 

Here = l/v is so called traversal time, which roughly corresponds to the time needed 
for a tunneling particle to traverse the classically forbidden underbarrier region. It is quite 




FIG. 2: Relative number of particles {N{t)/NQ) as a function of time (in units 1/r) for different 
values of the initial number of particles Nq/Nqp and uoXq = 3: diamonds — Nq/Nqp = 0; crosses 
— Nq/Ngp = 0.004; circles — No/Nqp = 0.04; boxes — Nq/Ngp = 0.4 



clear that at any reasonably high potential barrier the tunneling rate is small enough so that 
TtrN/N <^ 1. The ratio XN/Vq cannot be large and in many realistic cases it is even small at 
any reasonable number of particles in the droplet. We may conclude that the inequality (fT!?|l 
is robust and holds at all reasonable parameters of the system, which justifies the adiabatic 
approximation applied in this paper. 

Figure |21 shows the time dependence of the number of particle N{t) in a droplet for 
different initial numbers A^o- These numbers are scaled with the quantity Nqp defined in 
(jS)), which can be thought of as a number of particles, at which the nonlinear contribution 
to the single particle energy E{N) becomes comparable with the linear part E{0). We see 
that the time decay of the droplet is purely exponential, if the initial number A^^o is small. 
It deviates from this exponential behavior and becomes essentially slower at large initial 
numbers of particles in the droplet. This deviation is stronger the larger Nq is, but at large 
time when the total number of particles N{t) becomes smaller the behavior of the decay 
curve tends to the exponential one. 

We may also find (directly from (jH))) the initial slope of the decay curve as a function of 
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FIG. 3: Initial rate of evaporation F, in units of Fq as a function of the initial number of particles 
Nq/Ngp at for several choices of uoXq: diamonds — aoxo = 3 ; crosses — ooXq = 2 ; circles — 
aoXQ = 1 ; boxes — aoXo = 0.5 



the initial number of particles Nq, 

1 dN 



N 



NTn dt 



-aoXo 



1 + 2—— e " "^Vgp. (14) 
t=o V Ngp 

It is interesting to note that the nonlinear effects usually cause a decrease of the escape rate 
with the increasing initial number of particles Nq in the droplet as, e.g. is demonstrated 
in fig. El On one hand, this trend is determined by the exponential function in (fTlj) . 
corresponding to the density of particles at the exit point xq, and the decrease is stronger 
for higher and wider barriers, large axo- On the other hand the square root factor, i.e. the 
particle velocity at the exit point, may result in an increase of the escape rate. The latter 
may happen for rather low barriers when ao^^o < 1 (see fig- IS))- 

We have demonstrated here how the hydrodynamic approach to the description of a mul- 
tiparticle quantum system lead to a solution of the nonlinear problem of particle tunneling 
escape from a trap. Using the adiabatic approximation (neglect of Vt), justified for nonlin- 
earities, which are not necessarily small, we succeeded in obtaining an analytical solution in 
a one-dimensional case for the time dependence of the number of trapped particles and for 
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the tunneling rates at different nonlinearity values. The technique, we use, goes beyond the 
WKB approximation and, in particular, divergencies at the turning point of the trap poten- 
tial do not appear. That is why applying the same technique directly to three dimensional 
systems seems to be straightforward. Contrary to the standard WKB approach velocity is 
not necessarily a potential field (see, e.g. discussion in Ref. il6) and in three dimensional 
cases vortices can be considered. 

Acknowledgment Auhtors are indebted to S. Flach for discussions. 
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